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A Lattice Study of Spin and Flavour Symmetry in Heavy Quark Physics 
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We present a first test of heavy quark spin-flavour symmetries in matrix elements for semi-leptonic decays 
B and B DlV . We show that ©(l/mg) corrections are small at masses around the charm for form 

factors protected by Luke's Theorem, but are of order 30-40 % for hv and hv I^a-^ ■ 
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1. Introduction 

In the limit of infinitely heavy 6 and c quarks an 
exact spin-flavour symmetry relates the six form 
factors describing the semi-leptonic decays B —f 
Dlv and B —f D*li> to a single universal function, 
$,{v ■ v'), the Isgur-Wise (IW) function[l]. 

The form factors are defined as follows : 
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where v and v' are the four- velocities of the initial 
and final mesons. 

To be able to have any confidence in a determi- 
nation of the CKM-matrix element yci,[9][10] we 
need to reliably determine the leading symmetry 
breaking coefficients in an expansion in l/mq. 
The theoretical framework to include these cor- 
rections is provided by Heavy Quark Effective 
Theory. To parametrize the matrix elements at 
order l/mq one needs to introduce four addi- 
tional functions and a mass parameter[3]. To- 
gether with the IW function these functions de- 
cribe the interactions of the light degrees of free- 
dom in the heavy mesons and they are there- 
fore only calculable non-perturbatively. Typi- 
cally, Heavy Quark Symmetry (HQS) makes pre- 
dictions for the normalisation of decay rates or 
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ratios thereof at zero recoil, but the dependence 
on is not predicted. 

A lattice calculation of the meson weak decay 
form factors provides a systematic study of HQS 
from first principles. In this calculation we use 
the Wilson action with an 0(a) improvement [2]. 
A different approach can be taken by deriving an 
effective action where the heavy degrees of free- 
dom have been integrated out [4]. 

Recently, the complete expressions for meson 
weak decay amplitudes have been given to or- 
der l/mq and to next-to-leading order in QCD 
perturbation theory [5]. For the form factors of 
equation (1) we have 
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where (3i{mi^c,v ■ v') are short distance correc- 
tions. As we wish to display l/mq symmetry 
breaking effects we show plots of 
^ _ hi{v ■ v') 



1 + Pi{mi^c,v ■ v') 



(3) 



As defined in (2), (,riin{v ■ v') is renormalisation- 
group invariant and normalised to one at zero re- 
coil. 

2. Simulation details 

These preliminary results are obtained from 
60 quenched, /? = 6.2, 24^ x 48 lattices, cor- 
responding to a scale of = 2.73(5)GeV, as 
set by the string tension[6]. Quark propagators 
were generated using an 0(a)-improved action at 
two "bottom" hopping parameters, Kq = 0.129 
and 0.121, corresponding to spin averaged meson 
masses of 1.4 and 1.9GeV[7]; four "charm" hop- 
ping parameters, Kq, = 0.133,0.129,0.125,0.121, 
and three "strange" hopping parameters, k; = 
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0.14144,0.14226,0.14262. The vector(axial) cur- 
rents are 0(a)-improved versions of 67^(75)0. 
The method for calculating the matrix elements is 
well established [8]. The use of non-conserved lo- 
cal currents is taken into account by introducing 
rescaling factors Zy and Za- The determination 
of these constants is described in [9] and [10]. 

Errors are calculated using a bootstrap proce- 
dure on 100 samples and take into account corre- 
lations between different timeslices. 

By injecting momenta up to V2(;|^) at the cur- 
rent and momenta (0,0,0) and (1,0,0) for the B- 
meson we vary v ■ v' the range 1 to 1.3. 

3. Flavour symmetry - form factors at 
maximum recoil 

Our study of flavour symmetry is motivated by 
Luke 's Theorem which states that there are no 
terms of order l/mq in the hadronic matrix el- 
ements at zero recoil [3]. In terms of form fac- 
tors this means that only /i+(l) and hA^{i) are 
protected against first order corrections, i. e. : 
/Ij4j(1) = 1 + 0(^^). Corrections of order 

have been estimated to be less than 3% [5] and 
we may use these results to determine to good 
accuracy from either of these two form factors. 

To explore flavour symmetry we exploit the fact 
that knowledge of the IW function at maximum 
recoil, = 0, determines its dependence com- 
pletely. As we have computed three-point cor- 
relators over a wide range of b and c masses we 
may vary v ■ v' hy variation of the heavy meson 
masses[ll]. Define y by 

y.v.Aq^ = 0) = li^ + ^). (4) 
2 m£)» ruB 

For every mass combination we fit hA^{v.v') 
to Stech's relativistic-oscillator parametriza- 
tion(BSW) and interpolate to y. These values 
can then again be fitted to the BSW form which 
parametrizes the IW function in terms of the 
slope parameter —p^ = ^^(1): 

^,{vv') = ^^exp{-(2p2 - l)^i-4^}(5) 
1 + V ■ v' V ■ v' + 1 

We perform a two parameter fit to s(p[v-v') where 
the parameter s is introduced to absorb uncer- 
tainties in the overall normalisation. 

In Figure 1 we show a fit to the data for the 
heaviest of the light quarks, k; = 0.14144, a b 
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Figure 1. Isgur-Wise function, normalised to 
Z^^ . The slope is determined from a fit to data 
from different heavy quark masses at zero recoil. 
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Figure 2. Isgur-Wise function, normalised to 
Z^^ , at fixed heavy quark mass(circles) and from 
the data from Fig. l(crosses). 

quark of = 0.121 and all four values for the c 
quark. Fitting these points to s(p[v ■ v') we find 
= 1.4 1 4 whereas a fit to the data for the 
same light quark and degenerate heavy quarks of 
KQ = 0.121 yields p^ = 1.0 t 5(Figure 2). Find- 
ing good agreement we conclude that flavour sym- 
metry holds over the whole range of masses and 
that we are in the heavy quark limit at around the 
charm mass. Using the fact that Za is essentially 
mass independent[10] this allows us to determine 
p^ from a fit to the data for more than one mass 
combination. In Figure 3 we show a BSW fit to 
the data for the lighter of the two b quarks and 
all four c quark values. We find = 1.2 t I, 
consistent with the above results. 
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Figure 3. Isgur-Wise function, normalised to 
. The slope is obtained from a BSW fit to 
data from different heavy quark masses. 
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Figure 4. Ratio of hy/hA^ normalised to Za- 

4. Spin symmetry - ratios of form factors 

With the exception of and Ha^ the form fac- 
tors are not protected against 0(^^) corrections 
at zero recoil and the predictive power is very lim- 
ited. For hv we find deviations from the exact 
symmetry limit = 1 of order 30-40%, even 

though the lattice vector current is normalised 
using the relation (P, p|y^|P, p) = 2p^. 

Following Neubert[12], we study the ratio R = 
which can also be extracted experimentally. 
This ratio depends on one subleading IW function 
which, ultimately, we would like to extract from 
a lattice calculation. 

Here, we wish only to display the dependence 
of symmetry breaking effects for the ratio R at 
zero recoil. In Figure 4 we plot ZAR{i) for both 
b quark masses and all four c quark masses. We 



find large l/mq corrections to the exact symme- 
try limit P(l) = 1 of the order 30-40% (note that 
Za ~ 1.1 [10]). We see an increase of this ef- 
fect with decreasing b and c masses. This result 
is in good agreement with a sum rule prediction 
R{1) = 1.33±0.08[5]. 

Furthermore, we study the ratio h^/hA^ of the 
two form factors protected by Luke's Theorem. 
We find this ratio to be close to one at zero recoil, 
as expected, with a slight decrease for bigger v ■ 
v' . This is another indication that spin symmetry 
holds for masses around that of the charm quark. 

5. Conclusions 

We have shown that flavour symmetry holds for 
Ha^ for heavy quark masses around charm and 
that 0(^^) corrections to Ha^^ are small for v ■ 
v' < 1.3. 

The ratio h^/hA^ is close to one for v -v' < 1.3 
providing evidence for spin symmetry. 

Both hv and receive large 0(;^) correc- 
tions of order 30 to 40 %. 
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